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Abstract. - We study the quantum probability to survive in an open chaotic system in 
the framework of the van Vleck-Gutzwiller propagator and present the first such calculation 
that accounts for quantum interference effects. Specifically we calculate quantum deviations 
from the classical decay after the break time t* for both broken and preserved time-reversal 
symmetry. The source of these corrections is identified in interfering pairs of correlated classical 
trajectories. In our approach the quantized chaotic system is modelled by a quatum graph. 



A fundamental source of physical information are time-resolved decay measurements in 
open quantum mechanical systems. While the radioactive decay is a prominent paradigm, 
more recent experiments studied atoms in optically generated lattices and billiards [1-6], 
the ionization of molecular Rydberg states [7] and excitaton relaxation in semiconductor 
quantum dots and wires [8,9]. Most of these examples, and also the complementary theoretical 
investigations of quantum decay [10-20], address the semiclassical regime of systems with 
chaotic classical limit. 

However, despite this broad interest there is no satisfactory semiclassical theory for the 
observed quantum dynamics. It is known from numerical studies and random-matrix theory 
(RMT) calculations [16-19] that the quantum survival probability P{t) follows the exponential 
classical decay Pc\{t) only up to a break time t* . This break time scales with the number of 
open decay channels L and the classical life time tc\ as t* ~ ^fPtd [16]. For t > t* the 
quantum decay law is a universal function which depends only on L, tc\ and the Heisenberg 
time tH and is qualitatively different from Pc\{t) [17-19]. Up to now, none of these results 
was accessible by semiclassical calculations and thus their applicability to individual chaotic 
systems remained a matter of speculation. 

In this letter we show that a systematic semiclassical expansion for the quantum decay can 
be based on the van Vleck-Gutzwiller propagator. Specifically we obtain with this approach 
the above mentioned features for the quantum probability to survive inside a quantized net- 
work (quantum graph), which is one of the standard models in quantum chaos [21,22]. We 
identify the source of quantum deviations from the classical decay in the interference be- 
tween certain pairs of correlated classical trajectories (inset of Fig.P). We have calculated the 
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quantum corrections analytically to leading order in time for both, broken and preserved time- 
reversal symmetry (TRS). The resulting expressions, Eqs. H18|l . H19|) below, are in convincing 
agreement with the corresponding numerical data (Fig.^. 

Although we study in the present paper only systems with universal behaviour, the semi- 
classical approach which we develop bears the potential to include system-specific properties 
as well. It is natural that this important advantage over any ad hoc random-matrix assump- 
tion comes at the price that semiclas steal calculations cannot be completely independent of the 
underlying model. Therefore some technical details of our calculation are specific to quantized 
networks. Nevertheless our results are of interest also beyond this class of models because the 
same pairs of correlated trajectories will give rise to quantum corrections also in other chaotic 
systems. This expectation is based on the analogy to weak-localization corrections in spectral 
and transport properties [23-28] . 

It must be stressed, however, that the quantum corrections which we are describing here go 
beyond the known weak localization effects as their presence is not restricted to systems with 
time-reversal symmetry. In previous studies of pairs of interfering classical trajectories it was 
always found that they have no net effect on two-point correlation functions if time-reversal 
symmetry is broken. This applies to both, open and closed systems [23-28]. Corrections have 
been found in the case of shot noise [35], but this quantity is a higher-order correlator and 
involves groups of four correlated trajectories. So besides the new physical context in which we 
study the effect of classical action correlations the novelty in the present work lies in the fact 
that they can have a non- vanishing effect on two-point correlators even without time-reversal 
symmetry. 

The dynamics of a quantum system is governed by a time-evolution operator U propagating 
an initial state |^o) from time to t as 

IV't) - f/* IV'o) • (1) 

For a closed system U is unitary and therefore the norm of the initial state is preserved 
[['(/'tll^ = 1. In contrast, for an open system U is sub-unitary. In this case the norm of the 
state, i.e., the survival probability 

m? = ii'A^t) = {Mu'^u'\i^^) , (2) 

decays from unity to zero as a function of time. We describe the smooth part of this decay. 
Superimposed fluctuations are removed by considering an average P{t) = {\\ijjt\\'^)ipo,k over the 
initial state (V'o)- An additional average is over the energy k. Assuming that the dynamically 
relevant energy window contains G states (G — > oo in the semiclassical regime) and choosing 
a discrete basis |m), the average over It/jo) leads to 

P{t) = G-\tT [/*C/*t)^ = J2 ■ (3) 

m,n 

U can be approximated semiclassically by the van Vleck-Gutzwiller propagator [29,30] 

iU^)mn ^^Ap ejip{iRp{m,n;t)/h) , (4) 
p 

i.e., a matrix element describing the transition n ^ m has contributions from all classical 
trajectories p leading from m to n in time t. Rp is the action of the trajectory and Ap a 
complex amplitude combining a stability factor and an additional phase from the Maslov 
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Fig. 1 - The ratio between quantum and classical survival probability has been computed for a 
quantum graph with L = 10 attached decay channels and G = 500(+), 750(x), lOOO(n), 2000(o) 
internal states. It is shown as a function of the scaled time t = t/G. Quantum corrections are visible 
as deviations of the data from the horizontal line. For small t they follow (a) Eq. 118|l in the case 
with time-reversal symmetry and (b) Eq. 1191 without. These semiclassical predictions account for 
the interference within the pairs of classical trajectories shown schematically in the insets. The two 
trajectories forming a pair are identical along the segments a,b, . . . but differ in the crossing regions 
(solid vs. dashed arrows). 

index of p. With Eq. we obtain for the survival probability 

P{t) = ( ^ ApA; exp(i [Rp - Rg]/h)) , (5) 

where the summation is over all pairs of classical trajectories p, q which start and end at the 
same point and have not left the system up to time t. 

The classical probability to remain inside the system is reproduced if Eq. |(SJ) is restricted 
to the "diagonal" terms p — q. Then P{t) is simply a sum over all possible trajectories p with 
corresponding probabilities |-4pp. It decays exponentially as, 

P.iit) - exp(-t/tei) , (6) 

provided that the ergodic time of the chaotic flow is much shorter than the decay time td that 
is obtained from the relative phase-space area of the opening [31]. 

To justify the diagonal approximation, which was developed originally in the context of 
the spectral two-point correlator [32], one observes that in the semiclassical limit ?i — > 
the exponential in Eq. |(SJ) represents rapidly oscillating phases which cancel tipon averaging 
unless Rp^q are correlated. This is certainly the case for p = q, but recent work on action 
correlations [23,26,28,33-35] has shown that other pairs contribute as well. The correlations in 
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the actions stem from the fact that these trajectories are composed of long segments where they 
follow each other closely. In short crossing regions the order and/or orientation of the segments 
is modified to yield different trajectories p ^ q. It is convenient to represent trajectories by a 
symbolic code in which letters a,b, . . . stand for a whole segment. The codes of p and q are 
then different words composed of the same set of letters. The number of segments, i.e., the 
length of the code word determines the relative importance of the action correlations. 

Here we apply this scheme to the survival probability. As the trajectories p, q must start 
and end at the same points, their symbolic codes start and end with the same letter. We need 
then at least three segments to obtain different trajectories 

p = [abc] q = [abc] , (7) 

where b denotes the time- reversal of b (see inset of Fig. ^). These are precisely the pairs 
which are responsible for the weak localization contribution to the conductance [23] and they 
are present for TRS only. A contribution which is present for systems with and without TRS 
involves at least a four-letter word 

p — [abed] q — [acbd] , (8) 

(left inset of Fig. ^). However, it will be shown below that the contributions from all such 
orbit pairs cancel. Therefore, in order to capture the leading-order quantum corrections for 
systems without TRS, we need to consider trajectory pairs with five segments of the form 

p = [abode] q = [adebe] (9) 

(right inset of Fig. ^). All other permutations with a, e fixed are excluded since two segments 
can be combined into a single one such that a pair of the form results. 

The following calculation will be based on the trajectory pairs IZj)-®- It is plausible that 
they give the leading-order quantum corrections because the two partner trajectories deviate 
only in a minimum number of permutation points. This assertion is supported by our final 
result and also by related work on spectral correlations and transport [23,26,28,35]. 

The next step is to perform the summation in Eq. (jSJ. For this purpose we need explicit 
expressions for the amplitudes and the actions of the classical trajectories and therefore we 
consider a specific model system. We assume a quantized network as it is commonly used, 
e.g., in mesoscopic physics [36] and quantum chaos [21,26,37]. The discrete time-evolution 
operator on a network with B directed bonds has the form of a product of two B x B matrices 

U{k) ^ 5 X D{k). (10) 

Here Dmn{k) = Smn^^''''"^ is a diagonal matrix containing phase factors which describe the 
free propagation along the bonds of the network. Im {m = 1, . . . ,B) denotes the bond lengths 
which are chosen to be incommensurate in order to avoid non-generic degeneracies. The 
wavenumber k will be used for averaging, = limfc^oo /o'^ dfc' (•). In particular, we 
have 

(e"=('™-'"))fc=<5„„. (11) 

The matrix S fixes the topology of the underlying graph and the classical transition proba- 
bilities between its bonds, Pn^m = j'S'mnl^- Pn->m specifies a Markovian random walk on the 
graph which is the classical analogue of Eq. H1U|) . We make the simplifying assumption that 
all transitions have equal probability. 
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The phases of S are still to be determined. We consider first a closed system. Then S must 
be unitary. As a second condition we note that for preserved time-reversal symmetry the 
time-evolution operator is a symmetric matrix up to a unitary transformation [38]. This is 
satisfied if S" = S"^ ^ since then 13 2 JJ D~2 is symmetric. These two assumptions together with 
Eq. H12|l lead to a natural choice for the matrix elements, Smn — B^^^^ e'xjp{2Tnmn/ B) [37]. 
In order to break time-reversal invariance we have to destroy the symmetry of 5* but at the 
same time preserve the unitarity and Eq. 1)12(1 . This is achieved by a simple transformation 
exchanging neighboring rows of S, namely 5^^' = A5'(^^ with A,„„ — (5„i,„_(_i)". We have 
checked numerically that, according to spectral statistics, a closed graph with 5*^^'^) is a generic 
model for quantum chaos in the presence or absence of TRS, respectively (see also [26,37]). 

Finally we need to open the system. A standard way to do this is to restrict the time 
evolution to a subset Q oi G < B internal states by the projector n*^^^ with non-zero matrix 
elements Ilmn = 1 for m = n < G, 

f/ = n(^) J7n(^^ (13) 

Effectively, the remaining L = B — G bonds arc then perfectly absorbing and play the role of 
attached decay channels (leads) . Nevertheless the set of these bonds L influences the dynamics 
on C/, e.g., via the identity 

SmnS„^j^i = Snn' ~ 'S'mn5'„j„/ Vri, n (14) 

meg „gL 

expressing the unitarity SS'' = / of the closed system. 

The Heisenberg time of a network model is given by the number of bonds, — B [21]. 
As our final results scale with tn it will be convenient to represent them in terms of the scaled 
time T = t/B. Whenever r is used, we imply the semiclassical limit ijj — » oo. This limit is 
taken with r and L fixed and we will keep only the leading order terms in tn = B. 

Substituting Eq. I|10() into Eq. Q and expanding {U*)mn one finds that a trajectory p of 
length t on the graph is just a sequence of bonds po, . . . ,pt with pq ^ n and pt = m. Amplitude 
and phase are given by Ap — Sp^p^_^ . . . Sp^p„ and Rp/h = kilp^^ + . . . + Ip^) , respectively [21]. 
A considerable simplification results from the fact that due to Eq. (|11|) only orbit pairs with 
equal total lengths survive. Consequently the phase factor in Eq. (O is absent. 

For p = q we obtain P,,{t) - Y.p \A? = {G/Bf = (1 - L/B)\ i.e., Pci(r) = exp(-LT). 
This is equivalent to Eq. |(BJ with tc\ = B / L. For the summation over p we have used that 
according to Eq. ((12(1 the total probability of a trajectory of length t is -B~*, and that there 
are G*^^ such trajectories as each bond pi is summed over the whole set Q. 

In trajectory pairs composed of s segments, the total time t is the sum of the lengths of the 
segments and of a contribution from each of the s — 1 crossing points t = s — l-t-to+ifa + '- 
Besidcs this constraint, ta,tb, . . . can take any value > 1. This yields a combinatorial factor 
{IZi) (BtY^^ / (s — 1)!. Within the segments, amplitudes pair to classical probabilities, 
Eq. ((12(1 . as in the diagonal contribution above. Only the bonds right at the crossing points 
between segments must be treated separately as a phase difference between p and q occurs 
there. We denote the factor from all crossing points for the moment by and proceed 
first to the summation over the inner bonds of the segments. For a segment i of length ti, 
containing ti — 1 inner bonds, this yields G^^{G/BY' . The prefactor G~^ is absent for the 
first and the last segment as the first or last bond is not involved in phase factors and can be 
summed. Including the explicit factor G~^ form Eq. Q we have thus G^~'' {G / By~^^^ 
G^~'* exp(— Lr). Further, G^~'^ cancels B'^~~^ in the combinatorial factor above to leading 
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order in B. All terms together finally yield 

Pp,(r)=r^-V(s-l)!*P,^ci(r), (15) 
where the phases remain to be calculated. We have 

^acbd — X/ ^l3aS*g^S^fjSp^Ss^Sgi^ (16) 
^adcbe = ^2 ^ PaSg^S^^S*gSs^ S*j^S^gS*0 

where a, f3, f3, . . . stand for the last bond in a, the first in b, the last in 6, . . . and the summation 
for each index extends over the open graph Q up to some crucial restrictions which we discuss 
now. Consider a pair in the form Q with b = f3b' (3, assume P = P and define a' = a(3 and 
c' = pc. Then the very same orbit pair can also be represented as p = [a'b'c'], q = [a'b'c']. 
We exclude this ambiguity by imposing the restriction P p. Similarly, for the pairs ||SJ| we 
impose /? 7^ 7, /3 7^ 7 and for the pairs ^ P S, (3 S. With the help of these restrictions 
and Eq. f(T^ we have 

Khc = E E s,,si^s^^s;p (17) 

= E I E ~ E I ^PS,S*pa^^^S*^p. 

In the second term with P = p the amplitudes S combine into and thus this term yields 

L^G/B^ which is negligible for B 00. In the first term we apply Eq. (|14|l to perform 

the (unrestricted) summation over P,p. We obtain LG^ / B^ L [B ^ 00) plus another 
negligible term. Thus from Eq. p5|l with s = 3 we obtain 

Fg^(r) = (L/2)r2F,i(r). (18) 

Applying repeatedly Eq. p4|) and estimating the order of the resulting terms for _B ^ 00 we 
find further 

Pacbdir) = 0, 

Padcbeir) = (LV24)t4Pci(t). (19) 

These expressions agree with the RMT [17, 18] to leading order in r. This can be attributed 
to the rapidly mixing dynamics of our model, see Eq. \i'2.\ . It is expected that higher orders 
are reproduced as well, if in addition to Eqs. 0-10 trajectory pairs with more segments are 
included. Their omission is also responsible for the remaining deviations between Eqs. (|18|l . 
p9|l and our numerical data (Fig.0. These deviations are larger for TRS, where additional 
pairs with s = 4, 5 exist, while further corrections independent of TRS have at least s > 6. 

In conclusion, we have obtained semiclassical expressions for the leading-order quantum 
corrections which determine the dynamics of the survival probability after the break time t* . 
As a corollary we are also able to calculate the break time itself within the semiclassical theory: 
Using the definition c = Pq(t*)/Pci(t*) (c is some constant threshold for relevant deviations 
from Pel) and substituting Pq = Pd -I- P^^^^ for TRS and Pq — Pd -I- Padcbe for broken TRS we 
get t* ~ Bj^fTj = VLtci which agrees with the break times evaluated numerically in [16]. 
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